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Synopsis 

Histone variant Htzl substitution for H2A plays important roles in diverse DNA transactions. Histone chaperones Chzl 
and Napl (nucleosome assembly protein 1) are important for the deposition Htzl into nucleosomes. In literatures, it 
was suggested that Chzl is a Htzl-H2B-specific chaperone, and it is relatively unstructured in solution but it becomes 
structured in complex with the Htzl-H2B histone dimer. Napl (nucleosome assembly protein 1) can bind (H3-H4) 2 
tetramers, H2A-H2B dimers and Htzl-H2B dimers. Napl can bind H2A-H2B dimer in the cytoplasm and shuttles 
the dimer into the nucleus. Moreover, Napl functions in nucleosome assembly by competitively interacting with non- 
nucleosomal histone-DNA. However, the exact roles of these chaperones in assembling Htzl-containing nucleosome 
remain largely unknown. In this paper, we revealed that Chzl does not show a physical interaction with chromatin. In 
contrast, Napl binds exactly at the genomic DNA that contains Htzl. Napl and Htzl show a preferential interaction 
with AG-rich DNA sequences. Deletion of chzl results in a significantly decreased binding of Htzl in chromatin, 
whereas deletion of napl dramatically increases the association of Htzl with chromatin. Furthermore, genome-wide 
nucleosome-mapping analysis revealed that nucleosome occupancy for Htzlp-bound genes decreases upon deleting 
htzl or chzl, suggesting that Htzl is required for nucleosome structure at the specific genome loci. All together, these 
results define the distinct roles for histone chaperones Chzl and Napl to regulate Htzl incorporation into chromatin. 
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INTRODUCTION 



The histone variant H2A.Z, or Htzl in Saccharomyces cerevisiae, 
replaces H2A in nucleosomes and has diverse roles in eukaryotic 
cells [1]. Histone H2A.Z regulates transcription and it is partially 
redundant with nucleosome remodelling complexes [1]. During 
transcription initiation, H2A.Z establishes open chromatin near 
TSSs (transcription start sites) Genome-wide studies have re- 
vealed that two H2A.Z-containing nucleosomes flank an NFR or 
NDR [nucleosome-free (depleted) region] near the TSSs [2-4] 
and Htzl proximity to the TSS influences gene expression [5]. 
H2A.Z is enriched in both centromeres and telomeres [6-8], and a 
lack of Htz 1 can repress telomere-proximal genes [9] and abolish 
cell-cycle progression in yeast [10]. The molecular mechanisms 



for these phenomena are not very clear. It is known that H2A.Z 
is lost at promoters and is enriched at the centromere when mi- 
tosis begins in mouse trophoblast stem cells [11]. Furthermore, 
in differentiating murine ESCs (embryonic stem cells), Htzl de- 
position increases chromatin accessibility [12]. H2A.Z also has 
immunology-related functions, showing a specific distribution 
pattern around the genome variation sites in human CD4 + T 
cells [13]. 

In S. cerevisiae, because of the importance of Htzl, the mo- 
lecular mechanism for incorporating Htzl into nucleosomes is 
of critical interest. Depositing and evicting Htzl from chromatin 
requires the coordination of the SWR1 remodelling complex, 
specific histone modifications and histone chaperones. Zhang 
et al. showed that Htzl deposition partially relies on Gcn5, an 
HAT (histone acetyltransferase), and Bdfl, an SWR1 complex 



Abbreviations: ChlP-chip, chromatin immunoprecipitation and tiling microarray chips; GO, gene ontology; HA, haemagglutinin; HAT, histone acetyltransferase; Napl, nucleosome 
assembly protein 1; NFR or NDR, nucleosome-free (depleted) region; PCC, Pearson correlation coefficient; Pol II, polymerase II; TSS, transcription start sites; TTS, transcription 
termination site. 
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member that binds acetylated histones [3]. In a human cell line, 
purified Htz 1 -interacting complexes contained the components of 
the SRCAP chromatin-remodelling and TIP60 HAT complexes 
[14]. Additionally, the acetylation of H3K56Ac (histone H3 at 
lysine 56) alters the substrate specificity of SWR1 in exchan- 
ging H2A with Htzl [15]. Htzl incorporation is catalysed by 
the SWR1 complex in an ATP-dependent manner [16]. Import- 
antly, both H2A-containing nucleosomes and Htzl-H2B dimers 
specifically stimulate the ATPase activity of SWR1 [17]. 

Chzl is an important component of SWR1 -dependent Htzl 
replacement. It is a Htzl-H2B-specific chaperone that delivers 
Htzl for H2A substitution [18]. It is relatively unstructured in 
solution, but it becomes structured in complex with the Htzl- 
H2B histone dimer [19]. Two short alpha-helices and an exten- 
ded linker roughly follow the Htzl shape (PDB ID: 2JSS). The 
chaperone-histone Chzl-Htzl-H2B complex is mainly driven 
by attractive electrostatic interactions [20]. Our previous work 
suggested that Chzl regulates H2B ubiquitination and subtelo- 
meric anti-silencing which is independent of Htzl [21]. Chzl can 
facilitate the disfavoring property of Sptl6 to Htzl -containing 
genes [22]. 

Napl (nucleosome assembly protein 1) is also an import- 
ant component for Htzl replacement of H2A. Napl can bind 
(H3-H4) 2 tetramers, H2A-H2B dimers and Htzl-H2B dimers 
[9,16,23]. Napl binds H2A-H2B dimers in the cytoplasm and 
the complex shuttles into the nucleus [24—26]. Moreover, Napl 
functions in nucleosome assembly by competitively interacting 
with non-nucleosomal histone-DNA. Additionally, unincorpor- 
ated Htzl-H2B dimers form a complex with Napl [26,27]. These 
multiple roles make Napl very important for normal cell activity. 
In yeast, deleting the napl gene significantly increases atyp- 
ical histone-DNA complexes [26]. Xue et al. suggested that 
Napl genetically interacts with Pol II (polymerase II) and reg- 
ulates H3K9ac [28]. Histone chaperones are known to escort 
histone variants through complicated mechanisms [29]. Specific- 
ally, Chzl does not appear to interact with Htzl in the cytoplasm, 
rather Napl interactions with Htzl-H2B dimers maintain a sol- 
uble pool of Htzl for transport to the nucleus where Chzl -Htzl 
interactions can occur [30]. 

Despite the advancements mentioned above, the roles of Chzl 
and Napl in the deposition and eviction of Htzl remain largely 
unknown. Two important open questions are: (1) do both Chzl 
and Napl escort Htzl from the cytoplasm to chromatin; and (2) 
what roles do Chzl and Napl have in depositing and assembling 
Htzl -containing nucleosomes. In this study, we used ChlP-chip 
(chromatin immunoprecipitation and tiling microarray chips) to 
determine the genome- wide binding loci of Htz 1 , Chz 1 and Nap 1 , 
in both wild-type and chaperone-deletion strains of S. cerevisiae. 
We combined these binding location data, with new nucleosome 
maps built by next generation sequencing. Our data demonstrate 
distinct roles for Chzl and Napl in incorporating Htzl into chro- 
matin. Chzl does not bind to genomic DNA, whereas Napl and 
Htzl competitively bind genomic DNA. Deleting the chzl gene 
decreased Htzl binding. In contrast, deleting the NAP1 gene in- 
creased Htzl binding. Interestingly, the acetylation level of Htzl 
is dramatically increased in chzl mutant. The nucleosome occu- 



pancy of Htzl -bound genes decreased in htzl or chzl mutants. 
Our results elucidated the functional roles of Napl and Chzl in 
depositing Htzl on chromatin. 



MATERIALS AND METHODS 

Yeast strains and growth conditions 

The yeast strains that were used in this study are indicated in 
Supplementary Table SI. The genomic integration of genes en- 
coding C-terminal HA (haemagglutinin) fusions was performed 
using a PCR-based, one-step method for gene modification [31]. 
Yeast transformations were performed using a lithium acet- 
ate/polyethylene glycol-based method [32] . The strains were veri- 
fied by PCR analysis of the tagged gene loci and Western blot 
analysis of the fusion proteins. The strains were cultured at 30 °C 
in YPD medium [1 % (w/v) yeast extract/2 % (w/v) peptone/2 % 
(w/v) glucose] unless otherwise indicated. 

Chromatin fractionation assay 

The chromatin fractionation assay was performed as [33] with 
minor modifications. Briefly, after spheroplasting, cells were 
washed with the washing buffer. The pellet of spheroplasts was 
then resuspended in EB buffer. The cells were incubated on ice 
for 5 min for lysis with vortexing every 1 min to make the total 
cell lysate. Then 100 fi\ of cell lysate was placed on the top of 
200 fil of EBX-S buffer. The mixture was spun at 12 000 rpm 
for 10 min at 4°C. There was a while chromatin pellet, a clear 
sucrose layer, and a yellow supernatant fraction on top. Ten mi- 
croliters of the upper yellowish liquid was removed and kept as 
the supernatant fraction. The white chromatin pellet was washed 
and resuspended in the EBX buffer, and the SDS sample buffer 
was added to the protein sample. 

Western blot 

Proteins samples were separated by SDS-PAGE and transferred 
to nitrocellulose membranes. The nitrocellulose membranes were 
blocked with 5% (w/v) non-fat dried skimmed milk powder 
in PBST solution. HA-tagged proteins were detected with a 
1:1000 dilution of mouse monoclonal anti-HA antibody. Myc- 
tagged proteins were detected with a 1:1000 dilution of mouse 
monoclonal anti-Myc antibody. Rabbit polyclonal antibodies 
were used to detect PGK1 (1:5000 dilution), Htzl (1:1000 di- 
lution), HtzlK14 acetylation (1: 1000 dilution) and Chzl (1:2000 
dilution). 

RNA preparation and gene expression profile 

Yeast cultures were grown at 30 °C to a density of 1 x 10 7 cells/ml. 
The cells were collected and immediately frozen in liquid nitro- 
gen. Total RNA was isolated by hot acid phenol extraction. The 
total RNA was treated with RNase-free DNase I and purified 
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with a Qiagen RNeasy kit. Microarray labelling and hybridiza- 
tion reactions were performed [34] . Two-colour microarrays were 
performed using Agilent whole-genome 5. cerevisiae arrays. All 
experiments were performed with duplicate experimental and 
technical replicates of each condition. Genes with a ^ 2-fold 
change in expression compared with the expression in the relev- 
ant wild-type strains were considered to be significantly affected. 

Chromatin immunoprecipitation 

For each chromatin immunoprecipitation experiment, yeast 
strains were grown in the YPD medium to an OD^o of 1.0. ChIP 
was performed as described previously [35]. In brief, proteins 
were cross-linked to their respective DNA-binding sites with 1 % 
(v/v) formaldehyde for 1 h at room temperature. Chromatin was 
then disrupted by glass bead lysis and sheared to an average 
size of 400 bp. Sheared chromatin lysates were incubated with 
antibody-conjugated magnetic beads overnight at 4 °C with rota- 
tion. Following incubation, crosslinks were reversed in both the 
ChIP and whole-cell lysate fractions, and samples were analysed 
by DNA tiling microarrays. 

ChlP-chip 

For genome- wide binding analysis of HA-Htzl and Napl- 
9xMyc, ChIP was performed as described above. Linkers were 
annealed to the ends of the ChIP and input [WCE (whole cells 
extract)] DNA samples, and DNA was then amplified by PCR. 
Amplified DNA from the IP and input samples were labelled 
using a ULS aRNA Fluorescent Labelling DNA Kit (Kreatech). 
Labelled DNA from the ChIP and input samples were hybridized 
to S. cerevisiae whole-genome tiling arrays (4x44 k; Agilent). 
Data were extracted using Agilent Feature Extraction software 
and analysed with Agilent ChIP Analytics software (Agilent). 

Nucleosome-positioning analysis 

Wild- type, htzl A, chzlA and naplA cells were grown in the 
YPD medium to an OD 60 o of 1.0. All cells were then treated 
with 1 % (v/v) formaldehyde for 20 min and incubated with 
125 mM glycine for 5 min. Cell permeabilization, micrococ- 
cal nuclease digestion, protein degradation and DNA purifica- 
tion steps were performed as previously described [33]. DNA 
samples were treated with RNase A and were separated in a 
2% (w/v) agarose gel to assess the nucleosomal content. Bands 
corresponding to mononucleosomal DNA were extracted using 
a Qiagen gel extraction kit (Qiagen). Mononucleosomal DNA 
libraries were prepared and sequenced using an Illumina Gen- 
ome Analyzer II (Illumina Inc.) according to the manufacturer's 
instructions. 



Coordinates of TSSs, TTSs (transcription 
termination sites) and REG sites 

The genomic coordinates of TSSs and TTSs were retrieved from 
the Tables function of the University of California, Santa Cruz 



(UCSC) Genome Browser (http://genome.ucsc.edu) [36]. The 
experimentally identified DNA regulatory regions, transcription 
factor-binding sites, and regulatory variants were also retrieved 
from UCSC [37]; these data are called REG sites. The transcrip- 
tion rate data of yeast genes were retrieved from a previous report 
[38] . The data of TATA-box containing and TATA-less genes were 
previously published [39]. 

Binding profiles of Htzl and Napl 

Binding profiles for the whole genome were generated using 
the binding log 2 ratio signal from ChlP-chip experiments. The 
binding of any genomic region that is covered by a probe is 
represented by the log 2 binding ratio of the probe; the binding 
of a region that is not covered by a probe is zero. An Htzl- 
bound gene has at least one Htzl -bound probe that is identified 
within the region from 500-bp upstream to 300-bp downstream 
of the TSS. Htzl-bound probes are identified by two criteria: the 
central probe's binding ratio is >2 and the P value is <0.05, and 
the P value of either of the two neighbouring probes is <0.25. 
The Napl -bound genes were identified similarly. The numbers of 
Htzl-bound and Napl-bound sites are shown in Supplementary 
Figure S 1 . 

Occupancy patterns of Htzl and Napl 

The occupancy pattern of Htzl (or Napl) around special sites 
(such as TSSs, REG sites and TTSs) is represented with average 
binding ratio (log 2 ) by aligning the binding profile for each gene 
at the special site. 

GO (gene ontology) analysis 

To analyse the dynamics of Htzl binding, 5419 yeast genes were 
grouped into 40 gene function clusters according to their GO an- 
notation using GO Slim Mapper (http://www.yeastgenome.org/ 
cgi-bin/GO/goSlimMapper.pl). Then, the number of Htzlp- 
bound genes was counted in each cluster in the wild-type and 
deletion mutant strains, and the P value of the enrichment 
was calculated with a hypergeometric distribution. To com- 
pare the binding dynamics, the PCC (Pearson correlation coef- 
ficients) of the average binding ratio (log 2 ) profile for each 
gene cluster were calculated in the wild-type and deletion 
mutant strains. Additionally, we directly performed a GO en- 
richment analysis for Htzl-bound genes using GO Term Finder 
(http://www.yeastgenome.org/cgi-bin/GO/goTermFinder.pl). 

Identification of motifs 

DNA sequences of 60 bp around the centre sites of the bound 
probes were extracted and submitted to MEME to generate the 
motifs using the default settings of MEME [40], 

Nucleosome occupancy data 

The raw sequencing reads were mapped on the S. cerevisiae 
strain S288c genome using Illumina CASAVA software with 
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default ELAND alignment parameters. Only the uniquely 
mapped reads with less than 1 bp mismatch were used in further 
analysis. The nucleosome occupancy was obtained as follows. 
First, the length of each read was extended 73 bp in the 3' direc- 
tion. Then, the Waston-strand reads and Crick-strand reads were 
oppositely shifted 73 bp. Finally, the occupancy of each genomic 
site was calculated as the sequencing depth. The final nucleosome 
occupation is represented as the fold of the occupancy relative to 
the average occupation. 



Identification of nucleosome dyad sites 

The dyad coordinates of nucleosomes were determined with 
a wavelet transformed-based algorithm. The nucleosome occu- 
pancy signal was first smoothed with Daubechies wavelet at level 
4. Then, the peaks that satisfied the following two criteria were 
identified: ( 1 ) the height is > 1 .2 (the occupation is 1 .2-fold higher 
than the average occupation); and (2) the FWHH (full-width at 
half-height) is not <73 bp. 



Nucleosome occupancy patterns 

The nucleosome occupancy patterns near the special sites (TSSs, 
REG sites and TTSs) were represented by the average nucle- 
osome occupancy profile, which was calculated by summing the 
occupancy signal at each genomic site and then dividing 
the summated signal by the gene number [41]. 



RESULTS 

Genome-wide localizations of the histone variant 
Htzl and its chaperones Chzl and Napl 

Chzl does not physically interact with chromatin 
To resolve the roles of Chzl and Napl in depositing Htzl 
onto chromatin, we first used a chromatin fractionation assay 
to separate cells into the soluble protein fraction and insoluble, 
chromatin-bound fraction. We found that Napl directly associ- 
ated with chromatin and showed a high affinity for chromatin. 
However, we did not detect a Chzl signal in chromatin pellets 
(Figure 1 A). Because Chzl is a Htzl-H2B-specific histone chap- 
erone [18], we infer that Chzl does not physically interact with 
genomic DNA for facilitating to deliver Htzl-H2B to chromatin. 

We also performed genome- wide ChlP-chip analysis of Htzl, 
Chzl and Napl. For each of the three proteins, we obtained the 
intensity values [binding ratio (log2> and P value] of more than 
40000 probes that covered the entire yeast genome and construc- 
ted binding profiles for three proteins (Supplementary Figure SI). 
Our results indicated that Chz 1 did not bind chromatin (log 2 bind- 
ing ratio<l) (Figures SI and IB), which was consistent with the 
chromatin fractionation test shown in Figure 1(A). In contrast, 
Htzl and Napl had relatively high binding ratios (Figures S 1 and 
IB) particularly around TSSs (Figure IB). 



Napl overlays with Htzl in binding to chromatin 
The binding profiles of Htzl and Napl were strongly correlated 
(r= 0.8726; P value < 10~ 200 ) (Figure 1C), suggesting that the 
occupancy of Napl significantly overlaps with that of Htzl. We 
defined a gene as Htzl -bound (or Napl -bound) if there was at 
least one Htzl -bound probe (or Napl -bound probe) in the re- 
gion from 500 bp upstream to 300 bp downstream of the TSS 
(see the Experimental procedures section). Based on this cri- 
teria, 638 of the 839 targets for Htzl were shared with Napl (P 
value = 4.4x 10~ 38 ) (Figure ID), suggesting a significant over- 
lap in chromatin binding between Htzl and Napl. However, we 
cannot exclude the possibility that the different binding of Htzl 
and Napl in this study may be due in part to the different tags 
used in the ChlP-chip experiments (e.g. 9xMyc for Napl and HA 
for Htzl). 

Napl is enriched at the nucleosomes near TSSs 
To understand the organizations of Htzl and Napl at promoters, 
we calculated the average binding ratio profiles for Htzl and 
Napl by aligning the binding profiles for 5419 yeast genes at 
TSSs (Figure IE). The profiles for both Htzl and Napl showed 
a low density at ~80bp upstream of TSSs, corresponding to 
the NFR, and peaks at the two immediate sides of the trough 
that decrease to the background level (i.e. 0) at about 600 bp 
downstream of TSSs (Figure IE). Thus, Htzl and Napl have 
similar distributions showing both factors enriched at promoter 
regions. 

Htzl occupancy overlaps with the nucleosome 
occupancy, and Napl binds the genomic DNA that 
wraps around Htzl 

The intimate association between nucleosomes and Htzl led us 
to explore the relative position of Htzl to nucleosome dyads. 
We determined the nucleosome positions by the next genera- 
tion sequencing and identified the position of nucleosome dyads 
using a wavelet-based algorithm. In the first (i.e. — 1) nucle- 
osome upstream of the TSS, Htzl is ~ + 24 bp upstream of the 
dyad position; in the first (i.e. +1) nucleosome downstream of 
the TSS, Htzl is ~ — 37 bp downstream of the dyad position 
(Figure IF). At the centre of the NFR, the Htzl profile shows a 
low trough (Figure IF). Moreover, the position of Napl is very 
similar to that of Htzl (Figure IF). This suggests that the Htzl 
occupancy overlaps with the nucleosome occupancy. Since our 
ChlP-chip experiments were performed with sonicated genomic 
DNA, we cannot infer that all chromatin associated Htzl or Napl 
are in nucleosomes. Especially, the previous study has showed 
that Napl binding may inhibit non-nucleosomal binding of H2A- 
H2B to genomic DNA [26]. In conclusion, our data demonstrated 
that histone chaperone Napl can associate with genomic DNA 
wrapped with Htzl. We also noticed that the binding loci of Htzl 
relative to nucleosome dyad sites were different in the — 1 and 
+ 1 nucleosomes (Wilcox signed rank test, P value = 1 .4 x 1 0 ~ 20 ) 
(Figure IF). This was probably caused by the different rotational 
conformation of the two nucleosomes (Figure IF). 
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Figure 1 Binding of Htzl, Chzl and Napl on the yeast chromosomes by Southern blot hybridization and ChlP-chip 
experiments 

(A) Chzl does not physically interact with the chromatin. (B) The binding ratio (log 2 scale) profiles of Htzl, Chzl and Napl 
in the region of 130-140 kb on chromosome 1. The binding ratio refers to the ratio of the hybridization values for Cy3 
[protein (Htzl or Napl)-bound DNA] versus Cy5 (genomic DNA). In the bottom panel, the arrows indicate the transcribed 
regions, the arrow direction indicates the direction of transcription, and the labels on the arrow indicate the gene names. 
Htzl and Napl show a high binding signal, especially near the transcription start locus, indicating their binding on the 
chromatin, while the Chzl signal is lower than log 2 (2). (C) Htzl binding is highly correlated with Napl binding [r = 0.8726, 
P value = 1.0xl0~ 201 (t test)]. (D) Venn diagram shows the numbers of Htzl-bound genes and Napl-bound genes (see 
the Experimental procedures); the hypergeometric distribution P value is indicated for the overlapped gene number. (E) 
Binding patterns for Htzl and Napl in the vicinity of TSSs. The average binding ratio profiles surrounding TSSs for 5419 
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Htzlp and Naplp are negatively correlated with 
transcription frequency 

The enrichment of Htz 1 near TSSs is closely correlated with gene 
expression levels [3,4,41]. We tested the relationship between 
the average binding ratios of Napl and Htzl across TSSs and 
mRNA frequency using previously reported mRNA frequency 
data [38]. Yeast genes were sorted into five groups according 
to their mRNA frequency. A negative correlation (r< — 0.93, P 
value < 0.01) between the mRNA frequency and the enrichment 
of Napl and Htzl surrounding TSSs was observed (Figure 1G). 
The highly expressed genes exhibited a relative paucity of Htzl 
and Napl around TSSs (between — 200 and + 200 bp). In fact, 
the lack of Htzl and Napl was accompanied by the disassembly 
of the corresponding nucleosomes (see below). We measured the 
distance between the two peaks flanking TSSs and discovered 
that the distance between the two TSS-vicinity peaks for both 
Htzl and Napl decreased with increasing transcription frequency 
(r ^ — 0.60) (Supplementary Figure S2). This suggests that Napl 
falls off the chromatin with the loss of Htzl -containing dimers in 
regions of highly transcribed genes. It is notable that we observed 
a negative relationship between Htzl and transcription frequency, 
as Htzl incorporation into nucleosomes is a preliminary step in 
opening chromatin for transcription when the Htzl -containing 
nucleosome is unstable [5,12]. In yeast, Htzl is present in pro- 
moters as the default state, whereas it is recruited upon activation 
in mammalian cells. In both cases, however, the presence of Htzl 
at the promoter promotes the recruitment of RNAPII [42,43]. 



Genome-wide Htzl distribution in chaperone- 
deletion mutant strains 

To systemically uncover the effects of the Chzl and Napl chap- 
erones on Htzl deposition in chromatin, we performed protein 
expression test for Htzl, Htzl acetylated on K14, and H2A in 
wild-type, chzl A, napl A and chzl Anapl A, deletion strains. We 
also determined the genome- wide occupancy loci of Htzl in the 
deletion strains. 



Deletions of chzl and napl have opposite effect on 
Htzl deposition in chromatin 

We first carried out a chromatin fractionation assay to compare 
the levels of Htzl on chromatin in wild- type and mutant strains 
(Figure 2A). Surprisingly, Htzl signals from chromatin were 
dramatically reduced in the chzl mutant. Interestingly, the napl 



deletion strain had higher Htzl levels in the chromatin pellet than 
the wild-type strain. Furthermore, the chzl-napl double mutant 
also had increased levels of Htzl on chromatin like the napl 
single mutant, which indicated that Napl may have a dominant 
role regulating Htzl incorporation. In contrast, the levels of H2A 
in chromatin were unaffected in mutant strains. 

The result of the chromatin fractionation assay was further 
supported by ChlP-chip data (Figures 2B, 2C and Supplement- 
ary Figure S3). With respect to the Htzl binding ratio signal, a 
stronger correlation was observed between the napl A and chzl A 
napl A strains (Fisher transformation, P value = 0) (Figure 2B). 
In the chzl -deletion strain, only ~l/4 of Htzl -bound genes was 
the same as in the wild-type strain (Figure 2C). Additionally, 
89.6% of the Htzl-bound genes in the wild-type strain lost the 
Htzl in the chzl -deletion strain (Figure 2C), indicating that Chzl 
can distribute Htzl to specific genes. 

In contrast, the number of Htzl-bound genes in the napl- 
deletion strain was almost three times than that of in the wild-type 
strain (Figure 2C), suggesting that Napl represses the binding of 
Htzl to chromatin. Moreover, this effect was enhanced in the 
double deletion of chzl and napl (chzl A napl A) (Figure 2C). 
More than 63.3 % Htzl-bound genes in the wild-type strain were 
retained in the napl A and chzl A napl A strains. In short, deleting 
chzl not only decreased Htzl deposition in chromatin, but also 
changed the deposition loci, while deleting napl increased the 
Htzl deposition. 

To understand whether the association of Napl with genomic 
DNA is dependent on Htzl, we constructed Napl-myc htzl de- 
letion strains. Based on our ChlP-chip data, we chose eight loci 
with significantly binding with Htz 1 and Nap 1 . We performed the 
ChIP and RT-PCR to analyse the binding of Napl in wild-type 
and HTZ1 mutant (Supplementary Figure SI 1). Interestingly, we 
discovered that binding of Napl at these eight loci is dramatically 
decreased upon deletion of HTZ1, which clearly prove the notion 
that Napl -binding is Htzl -dependent. 



Htzl occupancy patterns near TSSs in chaperone 
deletion mutants 

We compared the average Htz 1 -binding ratio profile around TSSs 
in wild-type and chaperone deletion mutant strains. Regard- 
less of the strain, two peaks flanked a deep trough near the 
TSS (Figure 2D), showing similarity in the shape of the bind- 
ing patterns. About 400 bp from the TSS, the intensity of the 
average binding ratio profile for Htzl is sorted as 



genes of yeast are shown. The x-axis represents the distance to the TSSs; the y-axis represents the average probe binding 
ratio. Htzl exhibits a trough flanked by two peaks on either side; Napl binding is similar to that of Htzl, but it has a 
higher signal level, especially around TSSs. (F) Occupancies of both Htzl and Napl overlap nucleosome occupancies. 
The average binding ratio profiles relative to the - 1 and + 1 nucleosome dyad sites (left panel and right panel) and 
the centre of NFR (middle panel) are shown; the top panels indicate Htzl binding, and the bottom panels indicate Napl 
binding. (G) The average binding ratio profiles of Htzl (left) and Napl (right) for five gene groups with different transcription 
frequency (transcripts/hour). The transcription frequency data of yeast genes were retrieved from a previous report [26,53]. 
Correlation coefficient between the transcription frequency and the average binding ratio ( - 200 to 200 bp relative to TSS) 
is r= -0.93 (Fisher transformation, P value = 1. Ox 10~ 2 ) for Htzl and r= -0.97 (P value = 1.0x10 ~ 2 ' 9 ) for Napl. 
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Figure 2 Chzl and Napl have different effects on Htzl binding to chromatin 

(A) Binding of Htzl decreases upon deletion of chzl, while it increases after deleting either napl or both chzl and napl. 
Acetylation of Htzl at lysine 14 (HtzlK14ac) increases upon deletion of chzl. (B) Correlation of Htzl binding in deletion 
strains, the correlation coefficient (r) is indicated, and the P value of Fisher transformation is zero for the three coefficients. 
(C) Venn diagram of Htzl-bound genes in wild-type and gene deletion mutant strains. The Htzl-bound genes are identified 
by the criteria described in part of Experimental procedures. The hypergeometric distribution P-value is indicated. (D) The 
binding patterns of Htzl around TSSs of 5419 genes in the chzlA, naplA and chzlA naplA strains are shown. The 
y-axis represents the average binding log 2 ratio. The double-deletion mutant strain shows the highest average binding ratio 
level, and the chzlA strain shows the lowest average binding level. (E) Htzl is enriched at the - 1 and + 1 nucleosomes 
in deletion mutant strains. The Htzl-binding patterns are shown around the — 1 nucleosome dyad, the center of the 
NFR, and the +1 nucleosome dyad in the wild-type and gene deletion (chzlA, naplA and chzlA naplA) strains. The 
deviations of the Htzl binding peak to the nucleosome dyad axis are indicated in the table. 
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chzl A < wild-type < napl A < chzlA naplA, corresponding 
with the effects of Chzl and Napl on the binding of Htzl (Figures 
2A and 2C). This suggests that the chaperone deletions change 
the amount of Htzl in chromatin, but they have little effect on 
the Htzl binding pattern shape near TSSs. The genomic coordin- 
ates of REG sites (i.e. the DNA regulatory regions, transcription 
factor binding sites and regulatory variants) were retrieved [37]. 
We found that Htzl binding tends to flank REG sites, with the 
binding profile showing an increase from +800 bp and a clear 
valley starting at ~80 bp from the REG sites (Supplementary 
Figure S4). Htzl also shows a similar binding pattern near TSSs, 
which is consistent with previous studies [44,45]. 

We used the mutant strains to plot the Htzl -binding profiles 
around the — 1 and + 1 nucleosomes and the centre of NFRs. 
The Htzl occupancy overlaps with the nucleosome occupancy 
(Figure 2E), suggesting that most of Htzl located in the chro- 
matin. The binding loci of Htzl relative to nucleosome dyad sites 
are also different (P value ^ 4.0x 10 _23 ) between the — 1 and 
+ 1 nucleosomes (table on the right of Figure 2E). 



Binding dynamics of Htzl 

Dynamics of Htzlp binding and Htzlp-bound DNA 
sequence motifs 

Although the average binding ratio profile of Htzl near TSSs is 
similar in wild-type and deletion mutant strains, Htzl binding 
is decreased at the promoter and increased in the gene body 
after deleting the chaperones. We divided genomic DNA into 
three categories: promoter ( — 500 bp ~ TSS), gene body (TSS 
~ TTS) and other region. Deletions of the two chaperones caused 
an ~9 % decrease in number of Htzl -bound probes at promoters 
(Supplementary Figure S5) providing further evidence that both 
chaperones play roles in guiding Htzlp to the proper binding 
sites. 

To analyse the regulatory dynamics of Htzlp binding in the 
deletion mutant strains, 5419 yeast genes were grouped into 40 
gene clusters according to their GO annotation (http://www. 
yeastgenome.org/cgi-bin/GO/goSlimMapper.pl). The average 
Htzl -binding ratio profile of each cluster was calculated around 
TSSs. The PCC was calculated using binding profiles for the 
wild-type and mutant strains, showing the effects of gene dele- 
tions on the Htzl-bound genes for the 40 gene clusters. 



Ribosome genes had the lowest Htzlp-binding ratio levels 
and showed large changes in the binding profile across all dele- 
tion mutant strains (Figure 3A and Supplementary Figure S6A). 
Genes that were associated with nucleic acid binding TFs activity, 
DNA binding, phosphatase activity and enzyme activity had bind- 
ing profiles that were highly correlated between the wild-type and 
mutant strains (Figure 3A). Genes encoding nucleic acid-binding 
TFs activity showed a typical occupancy profile with Htzl en- 
riched at the two immediate sides of TSSs regardless of deletion 
implying that neither Chzlp nor Naplp affects Htzl regulation 
of these genes (Supplementary Figure S6B). 

Htzl intrinsically favors AG-rich DNA sequences 
We extracted the Htzl-bound genomic DNA sequences and cal- 
culated the binding motif of Htzl using the MEME tool. MEME 
generated three types of DNA motifs for Htzl and three for 
Napl. The primary Htzl motif was adenine- and guanine-rich 
(AG-rich) DNA sequences [Figure 3B (left column) and Sup- 
plementary Figure S7]. Importantly, the AG-rich binding feature 
was not disrupted in the deletion mutants, suggesting that Htzl 
intrinsically favours AG-rich DNA sequences. The primary Napl 
motif was very similar to that of Htzl, as would be expected by 
the findings that Htzl and Napl bind the same genomic loci in 
wild- type and supporting the hypothesis that Napl competitively 
binds Htzl-contianing dimers. 

We used GOMO to track the target genes of the motifs. Res- 
ults indicated that the target genes of the AG-rich motif associ- 
ates with transcription factor activity and transmembrane trans- 
port process [Figure 3B (left column)]. The secondary Htzl mo- 
tif is associated with ribosomal subunit constituents [Figure 3B 
(middle column)]. The second tertiary motifs are completely dis- 
rupted in the deletion mutants. This is more obvious in chzl 
mutants suggesting that Chzl partly contributes to the Htzl- 
binding specificity. 

Htzl favours TATA-less promoters [3], Our results showed that 
the average Htzl -binding ratio profile for TATA-less genes was 
higher than that for TATA box-containing genes in both wild-type 
and deletion mutants (Supplementary Figure S8). Because TATA 
box-containing genes are highly regulated and TATA-less genes 
perform more housekeeping functions [39], which generally have 
high expression levels, we inferred that Htz 1 enrichment at TATA- 
less promoters caused nucleosome instability and facilitated the 



Figure 3 Htzl preferentially binds purine-rich DNA sequences, and its binding can be changed by Chzl and Napl 

(A) Dynamics of Htzl binding in wild-type and deletion mutant strains. Yeast genes are classified in 40 clusters based 
on their functions using GO Slim Mapper (http://www.yeastgenome.org/cgi-bin/GO/goSlimMapper.pl). The heat map 
indicates the correlation coefficients of the average binding ratio (log 2 ) profile (from - 500 bp ~TTS) for each gene cluster 
between wild-type and deletion mutant strains. The bar plots show the P value ( - logio(P)) of the GO enrichment for each 
function cluster; in the bracket on the bar, the first number indicates the number of Htzl-bound genes of each cluster, 
and the second number is the percent of the Htzl-bound genes of the cluster relative to the total Htzl-bound genes. The 
identification of Htzl-bound genes is described in Figure 1(D). (B) Motifs of Htzl-bound and Napl-bound DNA sequences. 
The DNA sequences of Htzl-bound probes are entered into MEME to generate the binding motif using the default settings. 
The motif is then inputted into GOMO to find the putative target genes and to analyse their associated GO terms. The 
summed significant GO terms are listed below the motif if found; the number in the bracket indicates the number of 
significant terms using the significance threshold q value ^0.05. Abbreviations MF, BP and CC mean molecular function 
biological process and cellular component, respectively. 



© 2014 The Author(s) This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC-BY) (http://creativecommons.Org/licenses/by/3.0/) 567 
which permits unrestricted use, distribution and reproduction in any medium, provided the original work is properly cited. 



H. Liu and others 




Ht/lprncan binding counu Naplp mem biwhrtg count* 

C 3CC 1JX 3JK 



V, 1 


• •••9 3> ■-. 


Ate/ 


• • • * & ,i. : 


tadal 


• ••• 




asm 




• • • 1 3 ;■ 




•••••«•■■ 


Asw/ Ate/ 


• t • ® * ; ! 









»•••»«- v 




• • # § s. 




• .,- 


• • # • 9 * ' 


• A 


• ■* 


• • • \ 


• •••<$ •> - 


<■-. , 


• • • * $ * ■ 


• ••«$ ft * 






• • % 9 # * f - 











• I | 


• > 




• •• - 




■ 1 1 • t - 


• • 5* - 


• • S • 


• ••»»f - 






'•••*« ■■■ 


• 


0 




t • • # 


• •3 






t • • ft 






• ••»!£•;:■■ 





Wl 

Sete/ 
Acta/ Ate/ 

Ajmjji/ 
Anciprf Ate/ 

AsiiW 
Ajwr/ A/rrr/ 



Figure 4 Napl function is dominant in gene expression 

(A) Correlation heat map of gene expression between the htzlA, chzlA, naplA, htzlA chzlA, chzlA naplA and htzlA 
naplA strains. Gene expression is represented as the expression level fold change between mutant and wild-type strains. 
The correlation coefficient is calculated for 5419 genes. (B) A plot of gene expression versus Htzl mean binding. The Htzl 
mean binding is the average binding ratio in the promoter region ( - 500-0 bp relative to the TSS); the gene expression 
is the fold change between the mutant and wild-type strains. (C) Growth of yeast with deletions of chzl, napl, htzl and 
their combinations. 



disassembly of nucleosomes near TSSs during transcription ini- 
tiation, which is consistent with previous study [3], 



Transcriptome analysis 

Napl is dominant to Chzl in gene expression regulation 
Although the histone chaperones were proposed to be involved in 
the incorporation of histone variants into chromatin, we wanted to 
understand whether histone chaperones could regulate the tran- 
scriptome. We performed genome-wide transcriptome analysis 
of the htzl A, napl A, chzl A napl A, chzl A htzl A and napl A 
htzl A strains. The correlation analysis showed that the transcrip- 
tome profile of htzl was poorly correlated with that of both chzl 
and napl (Figure 4A). These data showed that histone chap- 
erones may regulate transcription and this role is independent 
of the association with the histone variant Htzl. Interestingly, 



Napl showed strong dominant roles in regulating transcription. 
The transcriptome profile of the napl chzl double-mutant strain 
strongly correlated with that of the napl single mutant, which 
supported the dominant role of Napl in regulating the incorpor- 
ation of the histone variant Htzl into chromatin. 



No significant connection between Htzl and Napl 
binding and the corresponding transcripts 
Because both htzl and napl deletion can lead to dramatic global 
changes in the transcription of many genes, we wanted to un- 
derstand whether there was a direct link between the binding 
of these proteins and the transcriptional changes that followed 
deletions. We observed that there was no statistically significant 
connection between Htzl and Napl binding and the correspond- 
ing transcripts (Figure 4B). 
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Figure 5 Deletion of chzl increases HtzlK14ac 

(A) Average binding ratio profiles of HtzlK14ac in the wild-type and chzlA strains. (B) A Venn diagram shows the number 
of acetylated Htzl-bound genes in the wild-type and chzlA strains. (C) Gene ontology analysis for acetylated Htzl-bound 
genes in the wild-type and chzlA strains. 



chzl and napl do not genetically interact, Htzl and 
Napl have partial functional redundancy 
We assayed genetic interactions between Htzl and histone chap- 
erones Chzl and Napl. As shown in Figure 4(C), chzl and napl 
did not show any genetic interaction. However, the napl htzl 
double mutant showed a growth defect upon a shift to 37 °C, 
suggesting some functional redundancy between Htzl and Napl. 
We also tested whether the genetic interaction between Htzl and 
histone chaperones will exist on different drug-containing plates. 
As showed in Figure 4(C), htzl single mutant was very sensitive 
to the presence of the DNA-damaging agents methyl mefhanes- 
ulfonate, hydroxyurea and formamide, which indicated that Htzl 
may have an important role in genome stability. However, the chzl 
and napl single mutants did not show a growth defect on these 
plates. Consistent with previous findings that napl genetically 
interacts with htzl, htzl napl double mutants grew even slower 
than the htzl single mutant on these drug-containing plates. As 
previously reported, functional role of Htzl is partially involved 
in the loss of Swrl activity. We also tested the sensitivity of 
the swrl mutant on these plates. Although the swrl mutant also 
showed sensitivity to DNA damaging agents, htzl mutant grew 
slower than the swrl mutant, which indicated that Htzl may have 
Swrl-independent functions to regulate genome stability. 



Acetylation of Htzl 

Deletion of chzl increases the acetylation of Htzl 
Interestingly, deleting chzl increased the acetylation level of Htz 1 
(ac-Htzl). Around the TSSs, the average binding profile of ac- 
Htzl in the chzl A mutant was obviously higher than that of the 
wild-type strain (Figure 5A). Furthermore, the number of ac- 
Htzl -bound genes of the chzl A mutant was also greater than 



that of the wild-type strain (Figure 5B). GO analysis indicated 
that histone binding-related and RNA binding-related genes were 
significantly associated with ac-Htz 1 when chzl was deleted (Fig- 
ure 5C). Perhaps this acetylation reflects a compensatory mech- 
anism to silence the effects of additional Htzl incorporation. 



Nucleosome-positioning dynamics 

Nucleosome occupancy of Htzl-bound genes decreases 
with htzl or chzl deletion but does not change with 
napl deletion. 

We determined the genome-wide nucleosome positions by next 
generation sequencing. The nucleosome occupancy near TSSs of 
5419 yeast genes is shown in Supplementary Figures S9 (A)-(D). 
The average read intensity near TSSs is shown in supplementary 
Figure S9(E). The typical distribution pattern was repeated in 
both the wild-type and mutant strains (Figure S9E). An NFR 
is shown around TSSs and REG sites (Figures S9E and S10). 
The NFR was flanked by two well-positioned nucleosomes and 
a phased organization of nucleosomes (nucleosome array) was 
observed downstream of the TSS. 

We identified 869 Htzl-bound genes in the wild-type strain. 
Although the gross levels of the average nucleosome occupancy 
near TSSs were not significantly different (Figure S9E), the nuc- 
leosome occupancy of Htzl-bound genes dramatically decreased 
with htzl deletion, especially for the —2, — 1, +1 and +2 
nucleosomes near TSSs (Figure 6A). In contrast, the nucleosome 
occupancy of non-Htzl -bound genes (Figure 6 A) did not signi- 
ficantly change in the wild-type and htzl A mutant strains except 
at the + 1 nucleosome (Figure 6A). Moreover, we also found 
that Htzl-bound genes had lower nucleosome occupancy in the 
chzl A mutant than in the wild-type strain (Figure 6B). These 
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Figure 6 The average nucleosome occupancy of Htzl-bound genes changes with gene deletion and transcription 
frequency 

(A-C) Nucleosome occupancy profiles for Htzl-bound genes in wild-type, htzlA, chzlA and naplA strains. The Htzl-bound 
genes are identified in the wild-type strain. (D) Nucleosome occupancy profiles for genes with different transcription 
frequency (transcripts/h). (E) Nucleosome occupancy profiles for genes with differential expression (expression>2-fold; 
expression < 1/2-fold). The fold change of expression is measured between the wild-type and htzlA strains. 



results indicate that H2A cannot compensate for Htzl deficiency 
to assemble nucleosomes at Htzl-bound genes. 

Unlike the chzlA or htzl A mutants, the nucleosome occu- 
pancy in the naplA mutant was not different from that in the 
wild-type strain (Figure 6C). This indicated that deleting htzl 
or chzl disassembled nucleosomes near the TSS of Htzl-bound 
genes, whereas deleting napl did not have an effect on the nuc- 
leosome occupancy of Htzl-bound genes. This can be explained 
by the effects of Chzl and Napl on Htzl deposition. There are 



Htzl -containing nucleosomes near the TSS of Htzl-bound genes. 
Deleting htzl results in the disassembly of the Htzl -containing 
nucleosomes, so we observed the decreased nucleosome occu- 
pancy. Deleting chzl can also decrease Htzl in chromatin, which 
decreases the nucleosome occupancy. Deleting napl increases 
Htzl deposition in chromatin. That is, Htzl in chromatin is suf- 
ficient to form Htzl -containing nucleosomes. Thus, the nucle- 
osome occupancy is not different between the wild-type and the 
napl A mutant strain. 
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A broader NFR is associated with highly transcribed 
genes 

We plotted the nucleosome occupancy profiles for the genes with 
different transcription frequency (Figure 6D). A broad NFR and 
a low nucleosome occupancy level are associated with highly 
transcribed genes in both the wild-type and mutant strains. This 
is consistent with previous reports [46,47]. We also noticed that 
the genes with a great change in expression (>2-fold or < 1/2- 
fold compared with wild-type strain) were associated with a broad 
NFR (Figure 6E). 



DISCUSSION 



Incorporating Htzl into chromatin is required for transcription, 
DNA replication, DNA repair and chromatin integrity [2,4]. 
The process involves a chaperone network [29]. Both Napl 
and Chzl are involved in the deposition of Htzl in chromatin 
[18,23,26]. We revealed the distinct roles of the two chaper- 
ones in escorting Htzl from the cytoplasm to the nucleus and to 
chromatin. 

First, according to previous reports [26,30] and our results, 
Napl is the Htzl chaperone both in the cytoplasm and in chro- 
matin: it is involved both in importing Htzl from the cyto- 
plasm to the nucleus and in incorporating Htzl into nucleosomes. 
Moreover, Napl controls nucleosome quality by eliminating non- 
nucleosome histone-DNA interactions [26,48]. Also, Napl is 
involved in both nucleosome assembly and disassembly [29,48- 
50]. In contrast, there is no evidence that Chzl is associated with 
Htzl in the cytoplasm [30]. Our finding suggested that Chzl does 
not directly interact with chromatin. Htzl-H2B dimers only in- 
teract with Chzl in the nucleus [18,30]. Secondly, our results 
suggested that Chzl facilitated Htzl deposition in chromatin al- 
though Chzl does not physically interact with genomic DNA. We 
observed increased Htzl levels in chromatin in the napl A strain, 
suggesting thatNapl represses Htzl deposition. This is partly due 
to two mechanisms: Napl can competitively bind to DNA with 
non-nucleosome histone-DNA interactions [26,48], and Napl 
can disassemble nucleosomes [29,48-50]. Napl is dominant to 
Chzl in Htzl deposition, chzl Anapl A having nearly identical 
phenotypes with napl A. We also discovered that Napl and Chzl 
appear to have little impact on nucleosome positioning because 
we only observed changes in nucleosome occupancy and not in 
nucleosome position. Nucleosome dynamics are probably me- 
diated by chromatin remodelling machines [49]. We found that 
the — 2, — 1, +1 and + 2 nucleosomes contain more Htzl than 
other nucleosomes. Moreover, the Htzl-containing nucleosomes 
disassemble if the Htzl supply is reduced (htzl A or chzl A). This 
suggests that Htzl deposition is specific for those Htzl -bound 
genes. That is, H2A cannot completely compensate for Htzl de- 
ficiency at Htzl -bound genes. Htzl -bound DNA sequences are 
adenine- and guanine-enriched. Finally, we observed that the as- 
sociation between the transcription frequency and nucleosome 
occupancy near TSSs was not changed in the deletion mutants. 



Deposition of Htzl on promoters has been showed to be typ- 
ical yeast genes whose transcription is frequently reprogrammed 
[51]. Two Htzl-containing nucleosomes flank the NFR, serving 
to poise quiescent genes for activation and transcriptional initi- 
ation [52]. No correlation between Htzl binding and changes in 
gene expression in HTZ1 mutant, and Htzl prefers to associate 
with low-frequency transcription genes further demonstrated the 
binding of Htzl appears to establish an epigenetic markers for the 
rapid activation of its associated promoters [35] . Thus, we thought 
that Htzl deposition at promoters is a preparation for transcrip- 
tion. The triggering of transcription needs other events, such as 
TF and Pol II binding. That is, alteration of Htzl binding will 
change the transcription preparation for quiescent genes. Sim- 
ilar to HTZ1 mutant, there is no significant correlation between 
Napl binding and deregulated transcription in NAP1 mutant. 
Histone chaperone Napl has been implicated to eliminate the 
non-nucleosomal histone DNA interaction [26]. It will be quite 
interesting to fully understand the functional significance of Nap 1 
binding. Additionally, as shown in Figure 3(A), kinase activity 
genes are enriched for Htzl binding when the chaperone was 
deleted. However, helicase activity genes and enzyme regulator 
activity genes were no longer enriched for Htzl binding in the 
deletion mutant strains, probably due to the overall increase in 
binding, whereas phosphatase activity genes were enriched only 
in the chzl deletion mutant strain, probably due to decreased 
overall binding. 

According to literatures and our present work, we proposed a 
model to illustrate the mechanisms of Napl and Chzl in depos- 
iting Htzl in yeast (Supplementary Figure S12). In this model, 
Napl and Chzl have different roles in Htzl deposition. Napl 
binds with Htzl-H2B dimer, forming a Napl-Htzl-H2B com- 
plex, in cytoplasm and then the complex is imported into the 
nucleus. In nucleus, Htzl-H2B-Napl interacts with SWR1 and 
completes the deposition of Htzl. Chzl does not associate with 
Htzl in cytoplasm. Htzl deposition mediated by Chzl is oc- 
curs in nucleus. During the nucleosome disassembly under con- 
ditions of transcription, DNA repair and other process, Htzl- 
H2B is released from chromatin into the nucleus. For the de- 
ciduous Htzl-H2B, their deposition need help of Chzl. Chzl 
catches the released Htzl-H2B and delivers it to SWRI and 
realize the deposition. However, there is another possibility to 
explain the increased level of Htzl incorporation in the absence 
of Napl. Rather than Napl and Htzl competing for binding 
to the same genomic loci, it may be likely that the absence of 
Napl may also decrease the incorporation of conical H2A, which 
would allow Htzl/H2B dimmers to assemble in inappropriate 
places. 
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